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(— I ' Context. This work examines the relation between the fractal properties of the photospheric magnetic patterns and those of the coronal 

magnetic fields in solar active regions. 
I ■ Aims. We investigate whether there is any correlation between the fractal dimensions of the photospheric structures and the magnetic 

O ' discontinuities formed in the corona. 

Methods. To investigate the connection between the photospheric and coronal complexity, we used a nonlinear force-free extrapola- 
C/3 , tion method that reconstructs the 3rf magnetic fields using 2d observed vector magnetograms as boundary conditions. We then located 

the magnetic discontinuities, which are considered as spatial proxies of reconnection-related instabilities. These discontinuities form 
well-defined volumes, called here unstable volumes. We calculated the fractal dimensions of these unstable volumes and compared 
them to the fractal dimensions of the boundary vector magnetograms. 

Results. Our results show no correlation between the fractal dimensions of the observed 2d photospheric structures and the extrapo- 
lated unstable volumes in the corona, when nonlinear force-free extrapolation is used. This result is independent of efforts to (1) bring 
the photospheric magnetic fields closer to a nonlinear force-free equilibrium and (2) omit the lower part of the modeled magnetic field 
volume that is almost completely filled by unstable volumes. A significant correlation between the fractal dimensions of the photo- 
spheric and coronal magnetic features is only observed at the zero level (lower limit) of approximation of a current-free (potential) 
magnetic field extrapolation. 

00 ' Conclusions. We conclude that the complicated transition from photospheric non-force-free fields to coronal force-free ones hampers 

' any direct correlation between the fractal dimensions of the 2d photospheric patterns and their 3d counterparts in the corona at the 

^\ , nonlinear force-free limit, which can be considered as a second level of approximation in this study. Correspondingly, in the zero and 

. first levels of approximation, namely, the potential and linear force-free extrapolation, respectively, we reveal a significant correlation 

between the fractal dimensions of the photospheric and coronal structures, which can be attributed to the lack of electric currents or 
to their purely field-aligned orientation. 
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1 . Introduction that the very nature of AR magnetic fields that indicates fractal 

analysis as a suitable tool for AR complexity determination. 
A variety of fractal analysis methods are available nowa- 

Solar active regions (ARs) have attracted the interest of many days. A basic classification of these methods distinguishes 

researchers over the years, since they are connected with solar between monofractal and multifractal techniques. An overview 

energetic events, such as solar flares and CMEs. The relation of of monofractal algorithms has been provided by McAteer 

AR magnetic complexity with their flare productivity has been et al. (2005), who indicates the perimeter area, the lin- 

lately widely investigated (Meunier 2004, McAteer et al. 2005, ear size area, and the box counting methods as the most 

Georgoulis 2008). There are many measures of the ARs' com- important monofractal methodologies. The pros and cons 

plexity, such as their size distribution (Harvey & Zwaan 1993, of these techniques, as well as the resulting discrepan- 

Meunier 1999), the length of the main polarity inversion line cies are discussed in the same work. These techniques 

(Falconer et al. 2006) or the magnetic flux along the polarity in- have been widely used in the literature: perimeter area 

version line (Schrijver 2007). A common approach nevertheless technique (Roudier & Muller 1987; Hirzberger et al. 1997 

relies on fractal analysis because the fractal dimension indicates Meunier 1999; Bovelet & Wiehr 2001; Janssen et al. 2003) 

the self-similarity of a structure over several size scales. The linear area technique (Tarbell et al. 1990; Lawrence 1991 

sizes of AR magnetic fields are known to display power-law dis- Schrijver et al. 1992; Balke et al. 1993; Meunier 1999) 

tributions (Harvey & Zwaan 1993, Abramenko 2005), meaning monofractal box-counting technique (Stark et al. 1997 
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Gallagher et al. 1998; Georgoulis et al. 2002), multifrac- 
tal box-counting technique (Lawrence & Schrijver 1993; 
Cadavid et al. 1994;Lawrence et al. 1996; Conlon et al. 2008) 
multifractal wavelet technique: (Hewett et al. 2008). This yields 
fractal dimensions for the 2d photospheric ARs {D2d) that vary 
from ~ 1.1 up to ~ 2.0, depending on the method used and the 
AR sample under consideration. 

The fractal dimension is also useful because it pro- 
vides a reliable test of various theoretical models and 
simulations against observations. Cellular automata mod- 
els (Isliker et al. 2000, Vlahos et al. 2002), percolation models 
(Seiden & Wentzel 1996, Fragos et al. 2004), random walk dif- 
fusion models (Lawrence & Schrijver 1993), and photospheric 
magnetoconvection simulations (Janssen et al. 2003) have been 
tested against observations through their fractal dimension. 

The majority of the above-mentioned fractal analyses has 
been implemented in 2d and applied to (mainly photospheric) 
magnetograms. Solar energetic events, however, take place 
above the photosphere. Because of the turbulent photospheric 
motions, the magnetic helicity, and shear generated in the coro- 
nal magnetic fields lead to magnetic discontinuities, which can 
give rise to nanoflares, microflares, flares, and CMEs when some 
critical threshold is exceeded. 

It would therefore be very interesting to seek the fractal di- 
mension Dm describing the complexity of the magnetic field 
in the 3<i volume above a given 7A magnetogram. Aschwanden 
and Aschwanden (2008a, 2008b) were among the first to fol- 
low this concept. They calculated the D2d fractal dimension of 
20 flares, using a standard box-counting technique on data from 
the Transition Region and Coronal Explorer (TRACE). To in- 
fer the DT,d in the corona, they created an analytical flare ge- 
ometry model. The coronal arcade geometry includes three free 
parameters (arcade length, width, and heUographic longitude), 
and makes a number of simplifying assumptions, such as (1) the 
arcade is near the equator and latitudinal projection eff'ects are 
neglected, (2) the magnetic shear along the neutral line is ne- 
glected, (3) the neutral fine is oriented in the east-west direction, 
(4) individual loops in the arcade are semicircular, and (5) the 
arcade length is assumed to be commensurable with the arcade 
width when observed near the limb. This work has shown that 
D2d and Dt,^ do not scale as Dxd/^. = £>3£//3, as expected for 
isotropic structures, but obey a more complex relation. 

The purpose of the present work is to elaborate the rela- 
tion between D2d and D^a further by using - for the first time 
- magnetic field measurements and extrapolations. We investi- 
gate whether a correlation between D2d and D^d actually exists 
and we have created an extended database of 38 ARs in the 
form of vector magnetograms. Our input magnetograms come 
both from the Imaging Vector Magnetograph (IVM; Mickey et 
al., 1996; LaBonte at al., 1999) of the Mees Solar Observatory 
and the spectropolarimeter (SP; Lites etal. 2001) of the solar 
optical telescope (SOT) onboard Hinode. We first calculated 
the D2(i of the 2d magnetic structures captured by the magne- 
tograms through a conventional box-counting algorithm. To in- 
fer the of the volume above a given AR, we extrapolated 
the observed photospheric magnetograms into the corona. The 
extrapolation method-of-choice was the nonlinear force-free op- 
timization method of Wiegelmann (2004), extending the work 
of Wheatland et al., (2000). Although nonlinear force-free ex- 
trapolation is still a fuUy open research topic, Wiegelmann's 
(2004) method has been recognized by comparison studies in- 
volving several nonlinear force-free extrapolation methods as 
one of the most reliable techniques (e.g. Schrijver et al. 2006, 
Metcalf et al. 2008). Potential extrapolation was also used as a 



zero level of approximation for reasons of comparison. After 
the extrapolation, we identified unstable volumes (UnVos) in the 
corona as ensembles of adjacent magnetic discontinuity sites. 
The 3d fractal dimension D^d is determined as the box-counting 
fractal dimension of these UnVos. 

This work is structured as follows. Section 2 describes the 
data used in this study along with the necessary corrections im- 
posed on them. Section 3 discusses the box counting and ex- 
trapolation techniques applied to our dataset, aiming to reveal 
correlations between 2d and 3d magnetic domains. Important 
parameters and quantities used are explained in detail, whereas 
work-around solutions on standard methodology drawbacks are 
suggested. Section 4 presents our results and discusses our find- 
ings. FinaUy, Sect. 5 summarizes our conclusions. 

2. Dataset 

NonUnear force-free extrapolation techniques typicaUy require 
vector magnetograms that are not as widely available as con- 
ventional line-of-sight magnetograms. Here we have created a 
database of 38 diff'erent AR vector magnetograms from the IVM 
and Hinode's SOT/SR 

The IVM obtains Stokes images in photospheric lines with 
7pm spectral resolution, larcsec spatial resolution (~ 0.55 arc- 
sec per pixel in full resolution) over a field of A.larcmirP' and po- 
larimetric precision of 0.1% (Mickey et al. 1996). We used both 
fully-inverted and quicklook IVM data. The quicklook data were 
obtained from the IVM Survey Data archive, made available at 
http : / /www. cora.nwra.com/ ~ ivm/ IVM -Survey Data /.The 
quicklook data reduction differs from the complete inversion in 
that it uses a simplified flat-fielding approach, takes no account 
of scattered or parasitic light, and no correction is attempted for 
seeing variations that occur during the data acquisition. On the 
other hand, the spectro polarimeter (SP) of SOT onboard Hinode 
obtains line profiles of two magnetically sensitive Fe lines at 
630.15 and 630.25nm, using a 0.16 X 164 arcsec sUt and has a 
spectral resolution of 30 mA (Shimizu 2004). 

In this study we used 10 fully inverted and 16 quicklook IVM 
vector magnetograms with 12 LevellD SOT magnetograms. To 
remove the intrinsic azimuthal ambiguity of 180° we used the 
Non-Potential magnetic Field Calculation (NPFC) method of 
Georgoulis (2005). For computational convenience we also re- 
binned the disambiguated magnetograms into a 128x128 grid. 

That our input data derive from instruments with diff'erent 
spatial resolutions should not influence our analysis, as the frac- 
tal dimension is scale-invariant and does not depend on the spa- 
tial resolution of the instrument used. This is true for mathemati- 
cal fractals however, because fractality holds for a finite range of 
sizes in reaUstic situations. For this reason, we first assume that 
all multi-instrument vector magnetograms belong to the same 
dataset, and then we duplicate the analysis using three different 
sets: one including fully-inverted IVM magnetograms, one in- 
cluding quicklook IVM data, and one including the Hinode/SOT 
data. 



3. Data analysis method 

Our analysis consists of 3 distinct steps. First we apply a box- 
counting algorithm to the 128 x 128 disambiguated photospheric 
magnetograms. We estimate the Dzd fractal dimension from the 
pixels with magnetic field strength exceeding a given thresh- 
old. Second, we apply the Wiegelmann optimization algorithm 
to our vector magnetograms in order to nonfinearly extrapolate 
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Fig. 1. Logarithmic plot of A'(e) versus e for AR 10953, where D^d = 
1.63. 

the magnetic field from the photospheric boundary. We thus con- 
struct a 3d 128 X 128 X 128 cube, within which the magnetic field 
is unambiguously determined. We also perform potential extrap- 
olation as a zero level of approximation. Third, we calculate D^j 
via box-counting the sites within our cubic grid which exceed a 
threshold of either the gradient of the magnetic field strength or 
the magnitude of the rotation (curl) of the magnetic field vector. 
The UnVos in the coronal volume are determined in this way. 
Finally, we investigate whether there is any correlation between 
the infeiTed D2d and D^d- The quality of the shown coiTelations 
is being judged by both the linear (Pearson) correlation coeffi- 
cients. 

3.1. Did Determination 

To calculate the Did of each AR on the photospheric level, we 
apply a standard box-counting technique, similar to McAteer et 
al. (2005). As input data we use the disambiguated 128 x 128 
vector magnetograms (see Sect. 2). The quantity examined is the 
magnetic field magnitude |B| at the photospheric level {z - 0), 
which can be unambiguously determined by the known magnetic 
field components Bx,By,B^ on the heliographic plane. A site 
is taken into account in the box-counting if its magnetic field 
magnitude exceeds a specific threshold B^r- If the number, A^, of 
strongly magnetized square boxes {\B\ > Bcr) scales with box 
size, e, as 

N( e) oc 6--°^" (1) 

then Did is the fractal dimension, whereby e is successively in- 
creased as e = 2", n = 0, 1, 2, ...7, with maximum emax - 128. In 
practice, Du is the scaling index of the power-law least-squares 
best fit between N{e) and e. The criterion \B\ > Bcr for identify- 
ing the photospheric ARs has been used widely in the literature 
(e.g. McAteer et al. 2005). 

However, the above-mentioned methodology is known to 
suffer from two major problems. First and foremost, it is sensi- 
tive to the selection of the critical threshold value Bcr- To deter- 
mine a non-arbitrary value of Bcr, we apply a histogram method, 
by constructing the histogram of the field magnitudes of all mag- 
netograms in our database. We then fit a Gaussian to this his- 
togram and define Bcr as the value, above which the histogram 



deviates from the Gaussian. This test yields Bcr = 230 G. The 
second drawback lies in the least-squares best fit being prone to 
large errors, mainly due to the small number of data-points in its 
dynamical range. To overcome this problem, we apply a good- 
ness of fit test in the form of a chi-square test, following Isliker 
(1992). The chi-square test is applied to a sliding window on the 
linear representation of the scaling relation, and it indicates with 
a 90% level of significance whether a range of the overall scaling 
is indeed a power law. In addition, we demand that the power- 
law scaling extends over at least one order of magnitude, in order 
to yield sufficient dynamical range for a reliable estimate of the 
fractal dimension. As an example in Fig.l we show the plot of 
log A^e versus log e for AR 10953, where Did - 1-63. 

3.2. Nonlinear force-free extrapolation 

The next step is to extrapolate the photospheric magnetic fields. 
Potential extrapolation is the lower limit (zero level of) approx- 
imation and it provides the simplest force-free magnetic config- 
uration. As such, it is expected to yield the best correlation be- 
tween Did and D^d- A linear, but non-potential, force-free field 
extrapolation is the first level of approximation, while the non- 
linear force-free (NLFF) field extrapolation is the second level 
of approximation. The most realistic treatment would be a non- 
force-free static extrapolation or a magnetohydrostatic / magne- 
tohydrodynamic model using the photospheric fields as bound- 
ary conditions. The latter, however, implies analysis and com- 
putational resources that far exceed the scope of this work. We 
chose the NLFF field extrapolation for two reasons: 

1 . While the linear force-free assumption may have some va- 
lidity in a minimum-energy AR corona, it cannot be trusted 
in the chromosphere or, even more, in the AR photosphere. 
There the relatively high value of the plasma beta-parameter 
implies Lorentz forces that cannot be neglected. 

2. Because the purpose of this work is to investigate whether a 
coiTelation between Du and D^d exists, it is crucial to refrain 
from using an extrapolation that is susceptible to artificial 
coiTelations between the 2d and the 3d domains. A linear 
extrapolation method would bind the 3d coronal magnetic 
field to its photospheric 2d boundary. 

The NLFF extrapolation used here is based on the op- 
timization technique introduced by Wheatland et al. (2000) 
and further developed by Wiegelmann (Wiegelmann 2004; 
Wiegelmann et al. 2006; Wiegelmann 2008). This technique re- 
constructs force-free magnetic fields from their boundary values, 
based on minimizing the Lorentz force and the divergence of the 
magnetic field vector in the extrapolation volume: 

L= r w(x,y,z)[\B\-\VxB)xBf- + \V-B\^]d\ (2) 
Jv 

In the above functional, w{x, y, z) is a weighting function and 
V denotes the extrapolation volume. A force-free state is reached 
when L = 0, w > 0. For w{x,y,z) - 1, the magnetic field must 
be available on all 6 boundaries of our cubic box for the op- 
timization algorithm to work. Nevertheless, real vector magne- 
tograms provide the magnetic field only for the bottom bound- 
ary, whereas the edge top and lateral magnetic field values re- 
main unknown. The weighting function is thus used to reduce 
the dependence of the interior solution on the unknown bound- 
aries. We introduced a buffer boundary of 100 grid points to- 
wards the lateral and top boundaries of the computational box. 
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Fig. 2. a) Potential extrapolation for AR 10953. b) The respective NLFF field solution, where the boundary vector magnetogram has been prepro- 
cessed. 



We then chose w(x, y,z) — I in the inner, physical domain and 
let VI' drop to with a cosine-profile in the buffer boundary to- 
wards the lateral and top boundaries of the computational box 
(see Wiegehnann 2004 for details). 

An additional useful attribute of Wiegelmann's NLFF field 
extrapolation code is the preprocessing alternative it offers. As 
the photospheric magnetic field is in principle inconsistent with 
the force-free approximation, a preprocessing procedure was de- 
veloped by Wiegelmann et al. (2006) in order to drive NLFF 
fields closer to a force-free equilibrium. Preprocessing mini- 
mizes the forces and torques in the system thus satisfying the 
force-free requirements more closely. While preprocessing the 
photospheric magnetograms, however, we applied some smooth- 
ing to the field components. To study the effects of this smooth- 
ing we duplicated the fractal-dimension calculation for unpre- 
processed photospheric fields, as well. As shown in Sect. 4, pre- 
processing - including smoothing - does not affect D2d signifi- 
cantly. 

Figure 2 depicts the field lines of the extrapolated magnetic 
field for AR 10953, as calculated by the optimization algorithm. 
Frame a) shows the current-free (potential) extrapolation for ref- 
erence while the respective NLFF field solution is shown in 
Frame b). 

3.3. Did Determination 

The fractal dimension D^d in the extrapolation volume is deter- 
mined by means of the box-counting and goodness-of-fit meth- 
ods described in Section 3.1. What differs in this case is the 
quantity used to determine which grid sites in the volume will 
be taken into account. This quantity wiU determine the UnVos 
present in the cubic grid. We have two different ways of calcu- 
lating these UnVos: 

1 . The average magnetic field gradient Gav 

For every site i, j, k within our grid, we calculate the average 
magnetic field gradient with its neighboring sites as 

^ _ \Bij,t-B,„^/nn\ 

^sum ~ i,^ sumxy^ sumyy^ sumz)' 

We define: 



Bv , , + By. ., , . 

-•i./.K+l <-l,]J<.-\. 



Depending on the location of each site within the volume, 
the number of nearest neighbors nn assumes the values nn = 
3, 4, 5, 6. The physical explanation for selecting this criterion 
lies in a steep gradient of the magnetic field strength being 
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thought to favor magnetic reconnection in 3<i, in the absence 
of null points (Priest et al. 2003). 
2. The normahzed magnetic field curl C„ 

For every site i, j, k within our grid we calculate the normal- 
ized magnetic field curl as 



This criterion will obviously highlight areas of high electric- 
current concentrations, which are known to play a role in the 
formation of magnetic instabihties. 

During the application of the box-counting method, we con- 
sider a site as unstable if it exceeds a critical threshold value (gra- 
dient or curl respectively). This critical value is determined with 
the histogram method described in Sect. 3.1, yielding Gcr - 0.09 
for the average gradient criterion and Qr = 0.005 for the nor- 
malized curl one. If the number. A', of unstable square boxes 
(Gov > Gcr or C„ > Ccr) scalcs with box size, e, as 



Nie) oc e-'^"' 



(3) 



then D-id is the 3d fractal dimension, whereby e is once again 
successively increased as e = 2", n = 0, 1, 2, ...7, with maximum 
Smax = 128. As an example in Fig. 3 we show the logarithmic plot 
of Nie) versus e for AR 10953, where D^d = 2.52 for the gradi- 
ent (frame a) and D^d = 2.91 for the curl (frame b). In both cases 
preprocessing were conducted prior to the NLFF extrapolation. 

4. Results 

As discussed in Sect. 3.2, several combinations were investi- 
gated with respect to the extrapolation methods and correlation 
coefficients. We first present the case of NLFF extrapolation with 
preprocessing. This is carried out before NLFF' extrapolation, as 
unpreprocessed magnetograms are not consistent with force-free 
extrapolations. Consequently it would be very hard to approxi- 
mate a nonlinear force-free equilibrium consistent with unpre- 
processed data. 

First we investigate the relation between D2d and Dj,d, when 
the latter is calculated through the average (Gav) magnetic field 
gradient. The corresponding Pearson correlation coefficient be- 
tween D2d and Djd is -0.154. The corresponding probabilities 
are shown in Table 1 for Pearson coefficient. We then use the 
normalized magnetic field curl (C„), to calculate D^d- In this sec- 
ond attempt the Pearson correlation coefficient between and 
Did is 0.251, as shown along with the corresponding probabil- 
ities in Table 1. Evidently, when NLFF extrapolation is used, 
there is no correlation between D2d and D^d at the 95% signif- 
icance level and this result is independent of the criterion used 
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Fig. 3. Logarithmic plot of N(e) versus 6 for the UnVos in AR 10953. UnVos have been calculated using Ga,, (Frame a) and C„ (Frame b). The 
respective fractal dimensions D^j are also shown. 



to quantify the UnVos. Moreover, independently of the magnetic 
complexity of an AR in the photosphere (as quantified by D2d), 
the complexity of the generated UnVos varies within a more or 
less well-defined range, as can be seen in Table 2. The absence 
of correlation between D2d and D^d, as well as the range of D^j 
values, is shown in frame a) of Fig. 4, where the NLFF case is 
presented. This figure suggests that - although the complexity in 
the photosphere may vary significantly - the UnVos in the corona 
retain a more or less well-defined behavior. 

If we attempt to locate these UnVos within the 128x 128x 128 
volume, one finds that more than 80% of them are accumu- 
lated in the lower layers, whereas only 20% are found at heights 
15 < z < 30, in units of the boundary magnetogram's lin- 
ear pixel size. The magnetic discontinuities are restricted to the 
lower coronal layers, whereas only a few weak additional dis- 
continuities are identified in higher z > 30 layers. This result is 
even better illustrated if we slice our volume in layers along the 
z axis and mark which sites per layer host magnetic instabilities. 
Figure 5 shows the results for AR 10953, where Go,, is used to 
determine UnVos in heights z = 17, 22, 27, 32 correspondingly. 
It is obvious that the higher they are from the photosphere, the 
less and the weaker they are. This can be shown alternatively 
by means of the fractal dimension: we identify which sites per 
layer Zo satisfy the relation Gav,j^,, ^ Gcr and by using the box 
counting technique we calculate the fractal dimension D^^ . The 

distribution of D^^^ with height is shown in Fig. 6. It is evident 
that each layer for < z < 20 is almost completely filled by 
UnVos, yielding a nearly Euclidian (non-fractal) dimension of 
^ 2. For heights 20 < z < 30, decreases from 2 to 1, 
while for z > 30 ZJj^ becomes smaller than 1, indicating very 
small and isolated, "dust-like", UnVos. 

This finding corroborates previous results. Regnier & Priest 
(2007) show that there is a clear preference for free energy accu- 
mulation very close to the photospheric level. Given this phys- 
ical property, we attempted an alternative approach by omitting 
the layers that are filled by UnVos and focusing only on the 
part of the volume where the filling with UnVos starts becom- 
ing sparse. In this case, would there be any correlation between 
the photospheric driver and the higher coronal structures? To ex- 
amine this possibility, we calculated D^d starting from the spe- 
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Fig. 6. Fractal dimension D'^ versus height z for AR 10953. 

Table 1. Pearson correlation coefficients between and D^j for the 
potential and the NLFF cases. 



Extrapolation 


Used 3d 


Correlation 


Corresponding 


Method 


criterion 


with Dij 


Probability 


Potential 


Gav 


0.719 


99% 


NLFF with preprocessing 


Gav 


-0.154 


65% 


NLFF with preprocessing 


c„ 


0.251 


87% 



cific z, above which only the remaining 20% of UnVos is identi- 
fied, thus excluding the lower coronal layers with a large filling 
factor. This leads to a decrease in the fractal dimension D^d, as 
expected, and an increase in the correlation coefficients, but we 
still find a lack of any significant correlation between D2d and 
Did- 

As the first row of Table 1 shows, the only case where we see 
a considerable correlation at the level of 95% significance be- 
tween D2d and D^d is when we use potential field extrapolation. 
Figure 4b illustrates the correlation between D2d and D^d when 
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Fig. 4. Fractal dimensions Dij and Djd for all 38 magnetograms in our sample. Frame a) depicts the case of NLFF extrapolation with preprocessing, 
whereas frame b) depicts the case of potential extrapolation. In both Frames a) and b) the solid line with squares represents in ascending order 
and the dotted line with circles stands for D^j when UnVos are identified by the G„y criterion, In Frame a) the solid line with triangles reflects D^^ 
when UnVos are identified based on the C„ criterion. 
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Fig. 5. Sites per z layer that host magnetic in- 
stabilities according to the Gav criterion for AR 
10953. a) z = 17, b) z = 22, c) z = 27, d) 
z = 32. 



Table 2. Range of D2J and D^j values, depending on the criterion used. 
100% of UnVos are considered. 



Extrapolation 


Fractal 


Used 


Range of 


Method 


Dimension 


criterion 


values 


Potential 


Did 


B„ 


1.57 + 0.43 


Potential 


Did 


G,„. 


2.331 ±0.331 


NLFF with preprocessing 


Did 


B„ 


1.64 ±0.34 


NLFF with preprocessing 


Did 




2.66 ±0.21 


NLFF with preprocessing 


Did 


c„ 


2.75 ±0.19 



the magnetic gradient criterion is used in the potential case. This 
is a reasonable result, considering that the potential extrapolation 
produces the simplest force-free magnetic configuration. At this 
zero level of approximation the significant correlation revealed 
between the photospheric and coronal structures is attributed to 



the lack of currents. Similar results can also be reproduced for 
the linear force-free extrapolation: the correlation between D2d 
and Did is again significant at the level of 99%, with a Pearson 
correlation coefficient of 0.75. Inspecting Fig. 4 and Table 2, 
we also notice that D^d is larger for NLFF field extrapolations 
compared to potential field extrapolations. This means more ex- 
tended UnVos, hence a higher spatial filling in the NLFF case, 
as expected. Other than that, UnVos show a strong preference 
for accumulating at low altitudes in the extrapolation volume for 
both NLFF and potential fields. Similar results were found by 
Vlahos & Georgoulis (2004), where linear force-free extrapola- 
tion is used. Their study does not apply a fractal analysis but 
identifies many magnetic discontinuities whose free magnetic 
energies and volumes obey well-formed power-law distribution 
functions. 
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Table 3. Pearson correlation coefficient of the fractal dimensions 
between the unpreprocessed and the preprocessed photospheric magne- 
tograms. 



Correlation 


Corresponding 


Coefficient 


Probability 


Pearson 0.863 


99% 



2.1n 




1 1 1 

5 10 15 20 25 30 35 40 
# Magnetograms 

Fig. 7. Fractal dimensions Dij for raw (dashed line with circles) and 
preprocessed (straight line with squares) magnetograms. 



Furthermore, we investigated whether the smoothing im- 
posed by the preprocessing has severely affected the magnetic 
field. For this purpose we compared the fractal dimensions D2d 
of the observed (unpreprocessed) and the preprocessed magne- 
tograms in the photospheric boundary. As shown in Table 3, the 
fractal dimension in the photospheric level is not significantly 
altered by preprocessing. The correlation between the prepro- 
cessed and raw magnetic data is significant at a level of 99%, 
as also illustrated in Fig. 7. While no evidence of correlation 
between D2d and 0^,^ exists in the NLFF field limit, significant 
correlation exists in the potential limit, even without preprocess- 
ing. 

Finally, we investigated how the spatial resolution differ- 
ences between the datasets derived from various instruments in- 
fluence our results. We separated the calculated fractal dimen- 
sions {D2d and D^) into three groups, depending on the instru- 
ment they are derived from. This separation substantially de- 
creases the number of data in each sample and is prompt to in- 
fluence our statistical results. Nevertheless, this is the only way 
to investigate how different spatial resolution and inversion so- 
phistication (in case of IVM data) can affect the analysis. Poor 
statistics will lead to relatively low confidence levels. To prop- 
erly define the threshold values for each of the subsets, we ap- 
plied the histogram test described in Sects. 3.1 and 3.2 for 2d 
and 3d correspondingly separately to each subset; nevertheless, 
the derived threshold values in all cases were very close to the 
global Bcr, Gcr, and Ccr values found for all 38 magnetograms, 
thus allowing us to retain them also for all separate subsets. 

Table 4 corresponds to Table 1 but only contains the fully 
inverted data derived from IVM (10 magnetograms). Table 5 
corresponds to Table 1 but only contains the data derived from 
HINODE (12 magnetograms). Table 6 coiTesponds to Table 1 



Table 4. Fully Inverted IVM subset: Pearson correlation coefficients 
between D2d and Dij for the potential and the NLFF cases. 



Extrapolation Used 3d Correlation Corresponding 

Method criterion with Dorf Probability 

Potential 0758 99% 

NLFF with preprocessing G„,, -0.278 61% 

NLFF with preprocessing C„ -0.102 25% 



Table 5. HINODE subset: Pearson correlation coefficients between D211 
and Dill for the potential and the NLFF cases. 



Extrapolation Used 3d Correlation Corresponding 

Method criterion with D^rf Probability 

Potential G^,. 0919 99% 

NLFF with preprocessing G„,, -0.5 91% 

NLFF with preprocessing C„ 0.427 85% 



Table 6. Quicklook IVM subset: Pearson correlation coefficients be- 
tween D2d and for the potential and the NLFF cases. 



Extrapolation Used 3d Correlation Corresponding 

Method criterion with D^rf Probability 

Potential 07765 93% 

NLFF with preprocessing G^v 0.238 65% 

NLFF with preprocessing C„ 0.099 30% 



Table 7. Pearson correlation coefficients of the fractal dimensions Dirf 
between the unpreprocessed and the preprocessed photospheric magne- 
tograms in the distinct data subsets. 





Data 


Correlation 


Corresponding 




Subset 


Coefficient 


Probability 


Fully Inverted IVM 


Pearson 


0.770 


99% 


HINODE 


Pearson 


0.989 


99% 


Quicklook IVM 


Pearson 


0688 


95% 



but only contains the quicklook data derived by IVM (16 mag- 
netograms). Finally, Table 7 coiTesponds to Table 3 but for each 
distinct dataset separately. Although substantial differences exist 
between the various datasets, it is evident that the results found 
for the total dataset are qualitatively retained also for the sub- 
sets of data coming from different instruments. The correlation 
between Dqj and D^j is significant for potential extrapolation, 
whereas it is absent in the NLFF limit. Preprocessing does not 
significantly alter the photospheric magnetic fields in any case. 

5. Discussion and conclusions 

This study investigates whether the complexity of the photo- 
spheric magnetic field correlates with the complexity of the 
UnVos formed in the corona. Using 38 vector magnetograms, 
we 

- estimate the fractal dimension D2d from the pixels with a 
magnetic field strength exceeding a given threshold Be- at the 
photospheric level through a standard box-counting method 

- extrapolate the magnetic field from the photospheric bound- 
ary using a 

1 . nonlinear force-free optimization algorithm with prepro- 
cessing at the photospheric level, 

2. standard potential extrapolation algorithm 

- calculate D^d by box-counting the sites within our cubic grid, 
which exceeds a threshold in 
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1. the averaged magnetic field gradient Gov with respect to 
their neighboring sites, 

2. the normahzed magnetic field curl C„ (only for the NLFF 
fields) 

- investigate whether there is any correlation between the in- 
ferred D2d and D^a- 

Our results show that there is no correlation revealed be- 
tween D2d and 0^,^ for the NLFF case. The spatial distribution 
of the UnVos with height shows that > 80% of the magnetic dis- 
continuities are accumulated in the lower corona (within 20 Mm 
from the photosphere). The system is evidently highly unstable 
at these heights, yielding processes that are clearly nonhnear. 
It is this strong nonlinearity at lower layers that does not allow 
the corona to respond proportionally to changes imposed by the 
photospheric driver or, conversely, does not allow the line-tied 
photosphere to strongly respond to a restructuring of the coronal 
magnetic fields caused by an energetic event (flare, CME). The 
results of Sudol & Harvey (2005) establish that, even for the 
largest flares, the photospheric response is significant, but rather 
small (median value of photospheric magnetic field change at 90 
G) even given that these flares have to occur close to the photo- 
sphere where most of the free magnetic energy resides. 

In contrast to the results for NLFF fields, we have obtained 
a significant correlation between D2d and D^d in the case of po- 
tential field extrapolations. As already explained, this should be 
attributed to the absence of currents, which forces the coronal 
magnetic structures to closely follow the photospheric driver. 
Therefore, the currents in the NLFF extrapolated fields lead to 
a more complex magnetic topology and the loss of correlations 
with the photospheric fields. 

It is interesting to notice the strong accumulation of UnVos 
close to the lower (photospheric) boundary in both the poten- 
tial and the NLFF fields. We believe this is not an artifact and 
that it has to do with the fine, fibril structure of the photo- 
spheric magnetic fields, which gradually fades as we move to- 
ward the corona. The structure of these forced fields is long 
known (Livingston & Harvey 1969; Howard & Stenflo 1972; 
Stenflo 1973) and it becomes less prominent at higher layers un- 
til the magnetic field fills the entire coronal volume, excluding 
small isolated areas of current sheets and tangential discontinu- 
ities in general (Parker 2004 and references therein). UnVos are 
designed to indicate these discontinuities and, as such, they will 
strongly accumulate at lower layers. This will happen regardless 
of the extrapolation method that will tend to create smooth fields 
in the volume, at the same time accommodating the finely struc- 
tured lower boundary. 

Further indications regarding of the absence of correlation 
between the photospheric and the coronal processes and struc- 
tures are provided by Metcalf et al. (1994), who examined the 
spatial and temporal relationship between coronal structures ob- 
served with the Soft X-ray Telescope (SXR) on board Yohkoh 
spacecraft and the vertical electric current density derived from 
photospheric vector magnetograms. Metcalf et al. (1994) found 
no evidence directly linking the electric currents observed in the 
photosphere to the heating of the coronal plasma indicated by 
the SXR brightness and temperature. 

The work of Aschwanden and Aschwanden (2008a, 2008b) 
- investigating the relationship between the fractal dimension 
of flare images captured by TRACE (in 2d) to the fractal di- 
mension of coronal arcades produced by an analytical geometric 
model (in 2>d) - indicates a complex relation between and 
DT,d. From our results, even this relation may be destroyed when 
fewer simplifications are used and the forced photospheric fields 



are taken into account. That the photospheric fields include sig- 
nificant Lorentz forces has been shown by Metcalf et al. (1995) 
and Georgouhs & LaBonte (2004). 

Our own and previous independent results seem to support 
the conjecture that the absence of correlation between the pho- 
tospheric and coronal fractal dimensions would stiU be the case 
(in fact, correlation should probably become worse) if a more 
realistic static or dynamic non-force-free modeling of the coro- 
nal field was used. Photospheric turbulence remains the driver 
for the coronal instabilities, but the strong nonhnearity of the 
system in the lower coronal layers destroys any kind of direct 
relation between the photospheric structures and their coronal 
counterparts. The photospheric driver forces the system to ac- 
cumulate a large number of magnetic discontinuities that store 
enough energy to explain the statistical properties of the solar 
activity in case of release. These discontinuities form patterns 
that do not follow the morphological properties of the photo- 
spheric magnetic flux concentrations, but have a strong impact 
on the expected dynamical activity of the system, namely, the 
magnetic energy release and the subsequent particle acceleration 
processes (Vlahos et al. 2004). 
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